Rotational structures of long-range diatomic molecules 
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We present a systematic understanding of the rotational structure of a long-range (vibrationally 
highly-excited) diatomic molecule. For example, we show that depending on a quantum defect, the 
least-bound vibrational state of a diatomic molecule with —C n /r n (n > 2) asymptotic interaction 
can have only 1, 2, and up to a maximum of n — 2 rotational states. A classification scheme 
of diatomic molecules is proposed, in which each class has a distinctive rotational structure and 
corresponds to different atom-atom scattering properties above the dissociation limit. 

PACS numbers: 33.15.Mt,34.10.+x,03.75.Nt,03.75.Ss 
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I. INTRODUCTION 

How fast can we rotate a molecule before breaking 
it |H 0? How does a rotational series terminates at 
the dissociation limit? How many rotational levels are 
there for a diatomic molecule in its last (least-bound), 
or next-to-last, vibrational state? These intriguing, and 
closely related, questions are taking on a new dimension 
of practical importance as our ability to make large sam- 
ples of long-range molecules (vibrational highly-excited 

molecules) ]3EI!ll,0ISH[l3, H3JE!_E3L_ arld even 

condensates of long-range molecules [lj, HH EH , contin- 
ues to grow. To understand the properties of a long-range 
molecule, especially how it responds to external pertur- 
bations such as collision with other atoms, we need not 
only the properties of a particular molecular state, such 
as the least-bound s state. We also need to know what 
are the states around it. It is this global structure of 
states that is the focus of this work. 

One approach to this problem is to compute the uni- 
versal spectra for each type of long-range interaction, 
—C n /r n , as we have done previously for n = 6 [IT], llq 
and n = 3 19], and simply observe them. This would, 
however, be very tedious and can never be completely 
inclusive. Our approach here is based on the recognition 
that the global structure of states, not including specific 
values for binding energies, depends only on the zero- 
energy wave function, more specifically, on its number 
of nodes as a function of both the angular momentum 
quantum number I and the exponent n characterizing 
the long-range interaction. 

Our results, and answers to the questions raised above, 
can be summarized in two simple formulas that will be 
derived later in the article. The first gives the dependence 
of the number of bound states on the angular momentum 
I for a quantum system with —C n /r n (n > 2) long-range 
interaction. 
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Here [x] means the greatest integer less than or equal to x. 
Ni is the number of bound states of angular momentum 
I. Nq is the number of s wave bound states. fi c is a 
quantum defect, to be defined later, that has a range of 
< /i c < 1. 

The second formula relates the quantum defect to the 
s wave scattering length 



a-Os 



a>0s 



tan(7r/i c ) + tan(-7r&) 
tan(7r/i c ) 



(2) 



Here ao s = ao//?n is the s wave scattering length, ao, 
scaled by the length scale j3 n — (2fj,C n /h ) 1 /(™ -2 ) asso- 
ciated with the long-range interaction; b = l/(n — 2); 
and 



a Qs — cos(irb) 



,2 b r(i-fe) 



r(i + 6) 



(3) 
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is the mean s wave scattering length of Gribakin and 
Flambaum |20f. scaled by (3 n . 

The consequences of these results are easily understood 
and are discussed in Sec. IIIII Equation QJ is derived in 
Sec. [n] It is another example of universal properties at 
length scale n , as discussed in a more general terms in 
two recent publications |2ll l22j . This universal property 
is followed by all molecules in varying degrees. Devia- 
tions from it and other issues are discussed in Sec. lIVI A 
primer of the angular-momentum-insensitive quantum- 
defect theory (AQDT) [171 Hs| , which is the foundation 
of this work, can be found in Appendix [S] 



II. DERIVATION OF EQUATION © 

Equation Q may be derived usin g tw o different meth- 
ods. One is to apply AQDT [13, llq . the version for 
arbitrary n > 2 as outline in [21j and Appendix^] to the 
zero-energy state of a diatomic molecule. This approach 
is discussed briefly in Appendix^ The other approach is 
the method of effective potential [2l|, I23 ■ It is this latter 
method that we present here, for the purpose of further 
promoting this powerful concept. While it makes no dif- 
ference in this particular case, for more complex systems, 



such as quantum few-body or quantum many-body sys- 
tems [22j where no analytic solutions are available, the 
method of effective potential may be the only way to un- 
cover universal properties at different length scales. The 
results would, of course, be mostly numerical in those 
cases. But a numerical solution done right can indeed 
yield universal behavior |22j . 

The method of effective potential is very simple, ft 
states that for a physical observable that depends only 
on states around the threshold, such as the number of 
nodes of the zero-energy wave function that we are look- 
ing at here, its universal behavior at length scale (3 n 
can be derived from any potential that has the right 
asymptotic behavior and is strongly repulsive at short 
distances. Specifically, a universal result at length scale 
/3„ is obtained from the corresponding result for the ef- 
fective potential by taking a proper limit to eliminate 
the shorter length scales while keeping the short-range K 
matrix, K c (0, /) fll7| and Appendix 1X1). to be a constant 
for one particular / J2J, |22| • 

We take here, for simplicity, a hard-sphere with an 
attractive tail (HST), 



^HST(r) 



c , r < r 

-C n /r n , r>r ' 



(4) 



as our effective potential. Its number of bound levels for 
angular momentum / is given by (|2l| and Appendix lA"}l 



JVhst(0 



m , j Uo ,m < 2/0 < ju ,m+l 
, 2/0 <ju ,l 



(5) 



where y Q = [2/(n - 2)](/3„/Vo) ( "~ 2)/2 , v Q = (2/ + l)/(n- 
2), and j Vo ,m {m > 1) is the m-th zero of the Bcsscl 
function J„ (x) |23| . 

Its K c parameter at zero energy is given by f|21| and 
Appendix [Ajl . 



*hst(0,0 



Jv {ya) cos(7ri/ /2) - Y Vo (j/o) sin(7T^ /2) 



Jv (2/o) sin(7ri/ /2) + Y Vo (y ) cos(7ri/ /2) ' 

(6) 
where J and Y are the Bessel functions 23] . 

In the limit of rn — > 0+ that eliminates the shorter 
length scale (see [23 for a more precise definition), yn> 
1, and the roots of the Bessel function are given by J2J] 

ju ,m^(m + Vo/2-l/4)7T. (7) 

K c (0,l) becomes an /-independent constant 

K c (0, l)^K c = tan(y + tt/4) . (8) 

Define the quantum defect, /i c (e, I), to be a parameter in 
a range of < /i c < f and related to K c by 



K c (e, I) = tan 



717/(6, Z) 



2(n - 2) 



(9) 



It is clear that fi c (0, /) also becomes an /-independent 
constant 
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where j is an integer chosen such that /i c falls in the range 
of 0</i c < 1. 

Combining these results, the number of bound states 
of angular momentum / can be written in the limit of 

r — > 0+ as 

JVhSt(0 -» TO . TO < J + M C - 732^ < m + 1 ■ (H) 
In other words, 

I 



JV HS t(0 
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where [x] again means the greatest integer less than or 
equal to x. Note that the result on the right-hand-side 
of this equation is now no longer just a property of the 
HST potential, but a universal property at length scale 
[3 n , applicable to any quantum system with the same 
long-range behavior and has a f3 n that is much longer 
than other length scales in the system. 

Since < /i c < 1, the integer j in Eq. I|I2|I is simply 
the number of bound states for / = 0. Eq. Q is thus de- 
rived. It is trivial to show that starting from a LJ(n,2n-2) 
effective potential [2l| (see also Sec. IIVB|) leads to iden- 
tical result. And again, the same result can be derived 
by applying AQDT to the zero-energy diatomic state, as 
outlined in Appendix lAl 

We emphasize that the method of effective potential is 
fundamentally different from the pseudopotential method 
[25j . In the latter method, a different pseudopotential 
is required for each partial wave. And each pseudopo- 
tential has at least one independent parameter to char- 
acterize the scattering of that particular partial wave. 
Without another theory relating scattering of different 
partial waves, no universal /-dependence of any kind can 
be established. (This is in addition to its well known lim- 
itations in describing bound states.) In contrast, a single 
effective potential is used to describe all /. And one can 
do this because a single parameter in AQDT describes a 
multitude of angular momentum states (see |l7l Il8j | and 
Appendix^). Put it in another way. Scattering of differ- 
ent angular momenta are indeed related, and so are the 
bound spectra of different angular momenta. Their re- 
lationships are determined by the long-range interaction 
([I3, 13 an d Appendix lA"|l . and are incorporated auto- 
matically in an effective potential |21|. The universal 
property described by Eq. (|TJ is but one reflection of this 
type of systematics. To further re-enforce this concept, 
we leave as an exercise to derive the relationship between 
the p wave scattering length and the s wave scattering 
length for n = 6. The derivation is trivial and all the in- 
formation required can be found in Appendix lAl and 26]. 
The answer will be provided in an upcoming publication. 



III. CLASSIFICATION OF MOLECULES USING 
QUANTUM DEFECT 

Let us first state that we do not include explicitly the 
effect of statistics when atoms are identical. It would 



TABLE I: Classification of diatomic molecules with —C n /r n 
(n > 2) long range interaction using quantum defect. Here 
Lmax.v is the maximum rotational quantum number for the 
vibrational state v. Av — v ma x — v. ao s = a,o/Pn is the scaled 
s wave scattering length. Note that some of the rotational 
states may be excluded for identical particles. Also note that 
scattering length has no definition for n = 3. 



Class 


Range of /j c 


-^max^v \^) 


Range of a,Q 3 





< fi c < b 


+ (n-2)Av 


2ao s < ao s < oo 


1 


b<fi c < 2b 


1 + (n - 2)A« 




71 — 3 


(n - 3)6 < fi c < 1 


n - 3 + (n - 2)Av 


— oo < aos < 



make our statements unnecessarily complex without in- 
troducing any new physics. In specific applications, all 
one needs to do is to exclude states that cannot satisfy 
the symmetry requirement (see, for example, |27j). as 
needed. 

The physical implications of Eq. Q) can be easily un- 
derstood by noting that No — 1 is the maximum vibra- 
tional quantum number, v max , while Ni — 1 is the maxi- 
mum vibrational quantum number, v max ^, that can sup- 
port a rotational state of angular momentum I. A vibra- 
tional state v can have all I for which v max .i > v. Letting 
L m ax,v to be the maximum rotational quantum number 
for vibrational state u, we have from Eq. Q 



1 



:L ri 



which can be written as 



L, 



(n - 2){v„ 



v) + [(n - 2V C 



(13) 



(14) 



where the square braket in the second term, [x], again 
denotes the greatest integer less than or equal to x 

This result suggests the classification of molecules into 
n — 2 classes, each corresponding to an equal interval of 
b = l/(n — 2) in the quantum-defect space. For class i 
with ib < /i c < (i + 1)6, we have i < (n — 2)/i c < i + 1, 
and therefore 



v,v(i) = i + {n- 2){v„ 



v) 



(15) 



Thus each class of molecules corresponds to a unique 
rotational structure that terminates at L max , v (i). This 
classification is summarized in Table [I] In particular, it 
means that the least-bound vibrational state can have 
1 (Class 0), 2 (Class 1), and up to a maximum of n — 2 
(Class n — 3) rotational states, depending on the quantum 
defect of the molecule. For the next-to-last vibrational 
state, add n — 2 rotational states to each class, and so on 
for lower vibrational states. 

What makes this classification useful is that each class 
not only has a distinctive rotational structure, it also cor- 
responds to distinctive atom-atom scattering properties 



above the dissociation limit. First, each class of molecules 
corresponds to a distinctive (non-overlapping) range of 
scattering length, which can be determined from Eq. @ 
and is summarized in Table [I] 

Equation J5J) derives easily from the definition of the 
mean scattering length 20], Eq. ©, the definition of 
the quantum defect, Eq. @, and the following rigorous 
relation between K c and the s wave scattering length 
(21] and Appendix 



«0.s 



>;2h 



T(l-b) 



r(i + 6) 



J ftT c (0,0) + tan( 7 r6/2) 
X c (0,0)-tan(7r6/2) 



(16) 



which is similar to the relation between scattering length 
and a semiclassical phase as derived by Gribakin and 
Flambaum [2{|. These equations combine to give 



n a.s 



a-o. 



tan[7r/i c (0,0)] +tan(7r&) 



(17) 



tan[7r/z c (0,0)] 

which is the exact relation between scattering length and 
quantum defect. It is the more rigorous way to write 
Eq. (0), applicable even when the system deviates from 
the universal behavior (see Sec. II V 01 and Appendix lAl . 

With the correspondence between quantum defect and 
scattering length, our classification of molecules can 
translate into other general statements, such as, a) The 
least-bound vibrational state of a diatomic molecule with 
o-Os > 2aos has only a single rotational state, b) The 
least-bound vibrational state of a diatomic molecule with 
negative scattering length has n — 2 rotational states. 
It is worth noting that molecules with negative scatter- 
ing length all fall into a single class, Class n — 3, while 
molecules with positive scattering length separate into 
n — 3 classes, from Class to Class n — 4. A similar 
feature was first noted by Gribakin and Flambaum }2(l ]. 

The different scattering properties for different classes 
are not restricted to the s wave. In fact, more interesting 
differences occur for higher partial waves. For example, 
Class does not have a p wave bound state for the last 
vibrational level. This p state, which would have been 
bound for /i c > b, does not disappear completely. It 
shows itself as p wave shape-resonance above the thresh- 
old, which actually becomes infinitely narrow (infinitely 
long-lived) as one approaches [f = b from the side of 
Class 0. In general, a Class i system is the one that has 
a shape-resonance of I = i + 1 closest to the threshold. 
The detailed properties of these resonances are however 
beyond the scope of this article (see, e.g. |17ll26ll2lll29|] ). 

The critical values of [i c = ib that are the boundaries 
between different classes correspond to having bound or 
quasibound states of angular momenta I = i + (n — 2)j 
(j being a non-negative integer) right at the threshold 
(Appendix 0. They have vibrational quantum numbers 
of v — Vmax — j, respectively This is a generalization of 
some of the results in [l|| to the case of arbitrary n > 3. 
Note that the wave functions for zero-energy bound or 
quasibound states are well defined and are given in the 
region of long-range potential by (Appendix FA)) 

(18) 



u e =oi(r) = Ar 1 / 2 J Uo (y) 



where r s = r/j3 n is a scaled radius, and y = [2/(n 



2)h 



-(n-2)/2 



This wave function has an asymptotic be- 



havior of 1/r at large r, thus representing a true, nor- 
malizable, bound state for Z > 0, and a quasibound (not 
normalizable) state for 1 = 0. The fact that s-wave wave 
function in the effective-range theory becomes completely 
meaningless when a^ = oo is only a limitation of the the- 
ory, not a reflection of any physical reality. 



IV. DISCUSSIONS 

We discuss here some special cases, deviations from the 
universal behavior, and how they might be treated. 



A. The case of n = 3 

Our result is applicable to n = 3, even though the 
scattering length has no definition in this case (for any 
Z) |30l l3l| . Specifically, it predicts that quantum sys- 
tems with n = 3 have only a single class (Class 0) with 
L m ax,v = v max — v. In other words, the last vibra- 
tional state for n — 3 has a single rotational state, an 
s state. The next-to-last vibrational state has two rota- 
tional states, and so forth. This prediction is confirmed 
by the analytic solution for —C^/r 3 potential (l9ll31| . 



B. The special case of LJ(n,2n-2) potentials 

For a set of Lennard-Jones potentials LJ(n,2n-2) (n > 
2) defined by 



V hJn {r) = -C n /r n + C 2n - 2 /r 



2n~2 



(19) 



the number of bound levels for any Z is given by ([2l| and 
Appendix ITj|l 



JV L Jn(0 



[z + \ - f ] , z > {u + l)/2 
, zq < (vq + l)/2 ' 



(20) 



where vq = (21 + l)/(n — 2), and the square braket [x] 
again means the greatest integer less than or equal to x. 
zo = (j3 n //3 2n ^ 2 ) n ~ 2 /[2{n-2)], where f3 2n -2 is the length 
scale assoicated with the C2n-2/r 2n ~ 2 interaction. Thus 
for the LJ(n,2n-2) potential, the universal dependence of 
the number of bound states on I, as specified by Eq. QJ, 
is exact, true even when f3 2n - 2 is comparable to f3 n and 
the corresponding potentials are so shallow as to support 
only a single or a few bound states. 

This result implies that to break the universal depen- 
dence on Z, one needs not only a short-range interaction, 
but the behavior of this interaction also has to be differ- 
ent from LJ(n,2n-2). 



C. Deviations from the universal behavior 

The key to understand qualitatively the deviations 
from the universal behavior is to recognize the origin 
of this universality. The universal Z-dependence origi- 
nates from the /-independence of A' c (0, Z) ([ItJ and Ap- 
pendix [SJ, which is a result of both the small mass ratio 
m e //i, where m e is the electron mass and /j, is the re- 
duced mass of the molecule (not to be confused with the 
quantum defect ^t c ), and the condition of (3 n 3> r$ where 
tq is a representative of other length scales in the system. 
(For HST potential, it coincides with the r^ that we used 
earlier.) 

With this understanding, it is clear that the universal 
behavior of Eq. QJ should be followed by all molecules to 
some degree. The mass ratio m e j ' \i is always small and 
can be taken for granted. (This is why we don't always 
mention it.) And almost by definition, f3 n is the longest 
length scale in the problem, otherwise it would not, and 
should not have been called the long-range interaction. 

It is also clear that the universal behavior is best 
followed by the states with highest vibrational quan- 
tum numbers. For example, consider our prediction of 
L max ,v(i) = i + (n — 2)(v max — v). For the least-bound 
vibrational state, v = v ma x, it would only require l- 
independence of K c over a range of AZ = n — 2. In 
comparison, the same result applied to v = v rnax — 9 
would require Z-independence of K c over 10 times that 
range, which generally becomes considerably worse (de- 
pending also on n, and other details of the short-range 
interaction). 

As far as predictions for the last few vibrational states 
(long-range molecules) are concerned, there is no need 
to worry about deviation except when \x c is very close 
to one of the critical values of /i c = ib, where a small 
Z-dependence may mean the difference between a bound 
state and a shape resonance. 

Whenever necessary, the deviation from universal be- 
havior can be easily handled within the AQDT frame- 
work. All we need is to count the nodes of the zero- 
energy wave functions more carefully! As discussed in 
Appendix^ AQDT is an exact formulation and an ex- 
cellent platform for exact numerical calculations. This 
also applies to node-counting: integrate the Schrodinger 
equation at zero energy and count the nodes up to a 
distance where K c (0, Z) has converged to a desired accu- 
racy [One computes K C (Q, I) by matching the integrated 
wave function to that given by Eq. I|A8(1 at different radii 
r. As a function this matching radius, K c (0,l) con- 
verges to a r-independent constant when the potential 
becomes —C n /r n and the wave function becomes that of 
Eq. (|A8|) ]. Adding to that the number of nodes beyond 
this distance, which can now be calculated analytically, 
gives one the total number of nodes. This way, one would 
never miss a node which could potentially be at infinity. 

One could also try to find if there are any systemat- 
ics in the deviation by going to the next, shorter, length 
scale. Any such attempt would however be necessarily 



system-specific and will be deferred to specific applica- 
tions. Examples of the universal rotational structure for 
n = 6 can already be found in [l3, uM > though they were 
not discussed explicitly. It was the simple structures ob- 
served there that motivated this work. 



V. CONCLUSION 

In conclusion, we have shown that the rotational struc- 
ture of a long-range molecule follows a simple univer- 
sal behavior that is characterized by two parameters, 
the exponent n of the long-range interaction —C n /r n , 
and a quantum defect, which is related in a simple way 
to the s wave scattering length [whenever it is well de- 
fined (n > 3)]. The resulting classification scheme gives 
a simple qualitative description of both the rotational 
structure of a long-range molecule and the correspond- 
ing atom-atom scattering properties above the dissocia- 
tion threshold. 

Finally, getting back to one of the questions at the 
beginning that we have not answered explicitly: how fast 
can we rotate a molecule before breaking it? The answer 
is, of course, L maXtV units of angular momenta, which is 
generally a very small number for long-range molecules. 
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APPENDIX A: AQDT: A PRIMER 

We give here a brief review of the angular-momentum- 
insensitive quantum defect theory (AQDT) |l7l Il8| . The 
focus will be on the conceptual aspects, and issues di- 
rectly related to this particular work. We point out that 
there are a number of different quantum-defect formu- 
lations for diatomic systems J22J |3a l^a |2J, 1^3 • There 
are also quantum-defect analysis |36l 1371 ] , and numerical 
methods that incorporate the concepts of quantum-defect 
theory 38]. Only our formulation is briefly reviewed here. 

Consider a radial Schrodinger equation 



2[i dr 2 



h 2 i(i + i) 

2fir 2 



V(r) 



u el {r)=0, (Al) 



where V(r) becomes —C n /r n beyond a distance r$. (tq 
does not have to have a precise value like the case of 
HST. It is a representative of length scales associated 
with interactions of shorter range, introduced to simplify 
our discussion.) 

In AQDT, the wave function in the region of long-range 
interaction (r > 7*0) is written as a linear superposition 



of a pair of reference functions 

u el (r)=A el [r eJ (r a )-K c (e,l)gl l (r s )} 



(A2) 



which also serves to define the short-range K matrix 
K c (e,l). The functions f c and q c are solutions for the 
long-range potential —C n /r n J3J, HJj- Their notations 
reflect the fact that with proper scaling and normaliza- 
tion, f c and g c depend on r only through a scaled radius 
r„ = r I ' [3 n and on energy only through a scaled energy 



7%_ =T/7 

Mm 



(h 2 /2n)(i/p n y 



(A3) 



Note that for the purpose of cleaner notion for arbitrary 
n, we have abandoned the factor of 16 used previously for 
n = 6 [13, EE, |2a |3!j , and the factor of 4 used previously 
for ra = 3 El HIl. 

Much of the art of a quantum defect theory [4Cll4lLl43 . 
I43| is in choosing f c and g c that best reflect the under- 
lying physics. For a molecule, the wave function at short 
distances is nearly independent of I because the rotational 
energy is small compared to electronic energy (originated 
from the small mass ratio m e /fi). AQDT takes advantage 
of this fact by picking a pair of solutions for —C n /r n po- 
tential that have not only energy-independent, but also 
/-independent behavior at short distances (possible be- 
cause n > 2) [HiJIil: 



9ti(rs) 



r</3„ 
r</3„ 



/ e c > s ) '-^" (2M^r^cos(y-7r/4) , (A4) 



-(2/7r) 1/2 r" /4 sm(2/-7r/4) , (A5) 



-(ra-2)/2 



where y = [2/(n — 2)] 

With this choice of reference pairs, matching to wave 
function at short distances yields an K c that is nearly 
independent of I, provided that tq is much smaller than 
P„ so that the reference functions at this point are well 
represented by their /-independent form of Eqs. (|A4(i - 
(IX5l) . 

An approximately Z-independent K c thus reflects the 
underlying physics that for long-range molecules and 
atom-atom scattering at low energies, the angular mo- 
mentum dependence is most important only at large dis- 
tances where its effects can be incorporated analytically. 
This is the physical origin of why a single parameter in 
AQDT is often capable of describing a multitude of an- 
gular momentum states [T^J UH |21j . The results of this 
work is but one reflection of the resulting systematics. 

The approximate energy-independence of K c under 
the same condition of f3 n ^> tq is fairly standard [4(1 142| . 
It is both because the reference functions have been cho- 
sen to be energy-independent at short distances, and be- 
cause the short-range wave function varies with energy 
on a scale of (7i 2 /2/i)(l/ro) 2 , much greater than the cor- 
responding energy scale associated with the long-range 
interaction, which is (?i 2 /2/j)(l//3„) 2 [l9| . 

In a multichannel theory that takes into account the 
hyperfine structures of atoms (starting from the formu- 
lation in J27| ), the concepts of AQDT, and the concept 



of /-independence in particular, remain unchanged and 
lead to an even greater reduction in the number of pa- 
rameters required for a complete characterization of the 
system |44ll4q 

We emphasize that AQDT is an exact formulation that 
does not require either the energy-independence or the l- 
independence of K c . It is simply the best framework, 
especially conceptually, to take advantage of them when 
they are there (/3„ 3> To). The parameterizations that we 
often use to extract universal behaviors should not dis- 
tract from the fact that AQDT also provides an excellent 
platform for exact numerical calculations, whether single 
channel J2l| or multichannel. This is especially true close 
to the dissociation limit, where matching to analytic solu- 
tions for ~C n /r n potential to obtain K c (e,l) converges 
much faster than matching to free-particle solutions to 
obtain the standard K matrix. The calculations in J21J 
are all based on this platform. [AQDT may even be the 
best method to calculate the scattering length when it is 
close to infinity] 

A major task of AQDT is, of course, finding the ref- 
erence functions. This is in general highly nontrivial 
[3l|,|33, especially analytically. No solution is yet avail- 
able for n = 5 at e ^ 0. This difficulty is however not a 
problem here as we need only the zero-energy reference 
functions, which can be easily found for arbitrary n and 



length. 



fl =ol (r s ) = [2/(n-2)] 1 /V s 1 / 2 [J, ( 2/ )cos(^ /2) 

-Y„ (y)sm(irv /2)}, (A6) 

9t=oi(r s ) = -[2(n-2)] 1 /2 r ]/2 [J ^ (y)sin(7ri/o/2) 

+Y l/0 (y) cos(ttv /2)], (A7) 

where v$ — (21 + l)/(n — 2). With these reference func- 
tions, the zero-energy wave function can be written either 

as 

u e =o*(r) = Mfl =ol (r s ) - K c (0J)gl =ol (r s )] , (A8) 

or as 

u e =oi(r) ^A' l r 1 J 2 [J l/0 (y)cos{a l ) + Y Uo {y)sm(a l )] , (A9) 

where ai = 7r[/i c (0,/) — lb] with the quantum defect 
/i c (e, /) being defined in terms of K c (e, I) by Eq. 10. 

The parameters K c and \x c both represent the same 
physics. K c is more convenient in computation, while 
fi c is able to represent all quantum systems in a finite 
parameter space of [0,1). In comparison, K c can take 
any value from -co to +oo. 

Equation @ is simply a result of node-counting the 
wave function, given exactly by Eq. (|A9() . from r to 
infinity (y = to j/o) 24J . Equation © is obtained 
simply by imposing the boundary condition u e =oi(r — 
r ) = 0. _ 

Equation l|16|) is obtained by comparing the asymptotic 
behavior the u e= oi(r) (for / = 0) at large r with the cor- 
responding expansion that defines the s wave scattering 



u e=0 i=o(r) —> A(r - a ) 



(A10) 



The derivation of Eq. JQ| in AQDT is straightforward. 
In the limit of rn <C /?„, the number of nodes of the 
zero-energy wave function inside ro is an /-independent 
constant [to a degree measured by the /-independence of 
K c (0,l)]. Counting the number of nodes of the outside 
wave function, Eq. I|A9|I . from tq to infinity (y = to 
Vo) HI- and ignoring the /-dependence of /i c (0, /) leads 
to Eq. Q. From this derivation, it is clear that deviation 
from the universal behavior is measured by the degree to 
which ^ c (0, /) or K c (0, I) is independent of L 

Having a bound or quasibound state right at the 
threshold corresponds to the boundary condition of 
u <L= oi(r) — > (a finite constant for / = 0) in the limit 
of r — > oo. Define 



xi(e) = tan[7r/i c (e, /) — nib] . 



(All) 



From Eq. I|A9() . the condition for a bound-state at the 
threshold is clearly 



xi(0) = tan[7r^ c (0, 1) - irlb] = 



(A12) 



which translates into ^ c (0,l) — ib for having bound or 
quasibound states of angular momenta / = i + (n — 2)j 
right at threshold, with corresponding wave functions 
given by Eq. JTSJ. In terms of K c , the same condition 
takes the form of 



K c (0,l) =tan 



n-2 



/ + 



1 



(A13) 



which is a generalization of the condition in |lj| to the 
case of arbitrary n. Note that the conditions expressed in 
the form of Eqs. (| A 12(1 and l|A13() are exact, with no as- 
sumption on the /-independence of either parameter. The 
universal behavior corresponds to when the /-dependence 
can be ignored (fi n ^> r ). 

The xi(e) parameter defined by Eq. (|A1 If) has also 
other applications. For example, for /^ c (0, /) >« ib, 
xi(0) >« is a convenient expansion parameter for de- 
scribing bound states of angular momenta / = i + (n — 2) j 
that are close to the threshold. xi(e) is also closely 
related to the Kf(e) matrix used in |26|. simply by 
Kf(e) = -xi(e). With this relation, all the results of [H 
can be rewritten in terms of cither K c or fi c . Making use 
of the /-independence of either parameter in these results 
leads, for example, to the relation between p and s wave 
scattering lengths. This has been left as an exercise, with 
the answer to be provided elsewhere. 

Reference [26| offers an important lesson on the im- 
portance of picking reference functions. The / and g° 
functions in [26|,|33, which define K°, differ from f c and 
g c only by trivial linear transformation. But because 
the resulting Kf(e) — — tan[-7r^ c (e, /) — irlb] did depend 
on /, relationships between scattering and bound spectra 
of different angular momenta were not recognized until 



much later [jj], [13 ■ Reference [2(| was only able to take 
advantage of the energy-independence of K® to show, 
for example, that the effective range and the scattering 
length are not independent, but are related in a way de- 
termined by the long-range interaction. The same con- 
clusion was also reached independently by Flambaum et 
al. using a different approach |47|. 



APPENDIX B: DERIVATION OF THE RESULTS 
FOR LJ(N,2N-2) POTENTIALS 

The analytic results for iV L j n (Z), Eq. IJ2H . and 
i4f£j n (0, 1) given in 21] [and therefore ^^(0,1), and of 
course ao s ] are derived from the zero-energy solution of 
the radial Schrodinger equation, Eq. IjAljI . for the class of 
potentials defined by Eq. (|19|l . Instead of giving all the 
boring math details, we prefer to simply note its relation- 
ship to the harmonic oscillator solution, as they have the 
same underlying mathematical structure. 



Upon a transformation x — (r/(3 n ) a and Ui(r) = 
x - ^ -1 ' ' ( 2a > Vi(x) with a = — (n — 2)/2, the correspond- 
ing equation at zero energy becomes 



tf#_ 



ft 2 7 ( 7 +l) 1 



+ —/ioj x 



E e 



vi(x) = 0, 



2/j, dx 2 2/ix 2 

(Bl) 
with 7+1/2 = [2/(n - 2)](/ + 1/2). Thus for 
the class of potentials given by Eq. (flfjfl - the solu- 
tion of the radial Schrodinger equation at zero en- 
ergy is equivalent to the solution of a 3-D isotropic 
harmonic oscillator with an effective angular momen- 
tum 7, a effective frequency determined by Tito — 
(fr 2 /2 M )(2/H)(/3 2 „_ 2 //3„)"- 2 (l//3„) 2 , at an effective en- 
ergy (not zero) E e — (h /2^)(l/a 2 )(l//3„) 2 . From this 
correspondence, both results are easily deduced. For ex- 
ample, the number of bound states is simply the number 
of harmonic oscillator states below, and including E e . 
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